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We use the dynamical mean-field method to determine the origin of the large ratio of the zero 
temperature gap to the transition temperature observed in most charge density wave materials. 
The method is useful because it allows an exact treatment of thermal fluctuations. We establish the 
relation of the dynamical mean-field results to conventional diagrammatics and thereby determine 
that in the physically relevant regime the origin of the large ratio is a strong inelastic scattering. 



Density wave ordering is a transition to a phase in 
which the electronic charge or spin density has lower 
symmetry than the underlying lattice. It occurs in a 
wide range of materials, including quasi-one-dimensional 
organic conductors , dicalogenides such as 2H — TaSe2 
|^],[|, 'A-15' materials like VaSi ||, 'blue bronzes' such as 
KM0O3 [Q, and the 'Verwey transition' material Fe3 04 
H . Much recent activity has related to actual or possible 
'stripe' density wave order in some members of the high 
temperature superconductor family Q and to charge and 
orbital order in some members of the 'colossal' magne- 
toresistance materials |^|. 

In most cases the density wave order evolves out of a 
metallic phase and in this situation it is usually believed 
[jlj to be driven by a fermi surface 'nesting' instability. 
The resulting equations are similar to those of the 'BCS' 
theory of superconductivity and in particular lead to a 
ratio of T = density wave gap A to density wave or- 
dering temperature T c which is close to the 'BCS' value 
Abcs/^bTc = 1.76. Almost all density wave materials, 
however, exhibit much larger A/k^T c ratios. In quasi- 
one-dimensional materials the large ratio may be under- 
stood as a consequence of critical fluctuations in low 
dimensionality, which decrease the transition tempera- 
ture more than the T — gap A. (Indeed, in a strictly 
one dimensional material, T c = while A > 0). However, 
A/T c values as large as 10 are also observed in quasi two 
dimensional systems and in many fully three dimensional 
materials [|],|||||| • These are not explicable either in the 
BCS approximation or in terms of the Migdal-Eliashberg- 
McMillan generalization || to so-called strong coupling 
superconductors such as Pb, which exhibit A/T c ratios 
only as large as 2.5. 

A generally accepted explanation in the physically rel- 
evant regime has not appeared. An extremely large 
electron-phonon interaction could localize carriers as po- 
larons with a large activation gap; the ordering temper- 
ature T c would then be controlled by a weak polaron- 
polaron repulsion and a large gap to T c ratio would result. 



However, this would imply a non-metallic normal state 
with a resistivity which is large and diverges rapidly as 
T — > 0, unlike the systems listed above (except perhaps 
for Fe304). McMillan made the intriguing suggestion 
(which has not been followed up by subsequent workers 
as far as we know) that the low T c (relative to A) was 
a phonon entropy effect jl0| : thermal fluctuations of the 
phonons would create areas where the local gap was small 
compared to the average, and a self consistent process of 
exciting the electrons into these areas would lead to the 
destruction of long ranged order. 

In this paper we address the issue quantitatively using 
the dynamical mean-field method | p2[ in the implemen- 
tation described in [ pH} . The technique permits a com- 
plete solution at all coupling strengths and temperatures, 
within a local (momentum independent self-energy) ap- 
proximation which is now generally accepted as reliable 
in d = 3 spatial dimensions jljj . We use it to show that 
in the weak to intermediate coupling regime relevant to 
most of the materials listed above, the crucial physics is 
inelastic scattering of electrons by the phonons involved 
in the lattice distortions. The phonon entropy effects 
proposed by McMillan |l(| are present, but are found be 
less important. This paper builds on the previous work of 
Ciuchi and DiPasquale who extended the methods of |l3} | 
to the density wave orderng case and computed the phase 
diagram in the intermediate to strong coupling regime 
Jli]] . We focus here on weak to intermediate couplings 
and present new physical interpretations. 

We study the Holstein model O], defined by 

ij i i 

(1) 

Here the operator c\ creates an electron on site i and 
Ti is the displacement of the ion on site i measured from 
the equilibrium displacement corresponding to to a uni- 
form distribution of electrons and rescaled so that the 
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electron-phonon coupling is unity. We have written the 
rescaled phonon stiffness in terms of an electron band- 
width parameter W defined below and a dimensionless 
coupling A. Because it contains all of the physics of rele- 
vance to us and simplifies the computation substantially 
we have adopted a classical limit (phonon kinetic term 
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Si 2M i s missing). We comment below on the effects of 
the neglected quantum fluctuations. We also assume a 
bipartite lattice, which for our purposes is a lattice pos- 
sessing a dispersion (fourier transform of ti—j) and 
a wavevector Q such that ££_q — f° r a U k in the 
Brillouin zone. An equivalent definition is that the lat- 
tice may be divided into two sublattices A and B, such 
that ti-j only connects A sites to B sites. We further 
specialize to n = 1/2; this implies a = 0. 

At n = 1/2 the ground state of (Q) may be shown to be 
charge ordered for all A by minimizing (if hoi ) over . To 
compute (ifhoi ) we introduce operators a) and creat- 
ing electrons on the A and B sublattices respectively and 
make the ansatz that = A on the A sublattice and = 
—A on the B sublattice. In momentum space the Hamil- 
tonian may be written as the 2x2 matrix ifhoi (A) = 
\N % + \ £5 (a* bQ [A r z + e % r x ] where r 2 and 

t x are the usual Pauli matrices and all k-sums are over 
the full zone. The resulting energy dispersion is Et = 

± y/e| + A 2 and the energy of Hhoi (A) may easily be 

found, and minimized over A. The equation yielding the 
minimum value is 
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If fj, = the integral, (||), is logarithmically divergent 
as tk — * and thus a A > solution exists. Later 
in this paper we use the semicircular density of states 
Pa(e) = i Efc^( e - H) 



-\JW 2 — e 2 for which 
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the equation can be solved analytically, yielding A = 
4|«: e -*(l + 0(A)). 

We now turn to the treatment of non-zero tempera- 
tures, following refs |T^,[l4|]. The crucial object in the 
analysis is the electron Greens function, which in the 
a — b basis introduced above is a matrix G given by 



G r 1 = iu n +fi— Y, n (iuj n ) - e fe r x - £ co (i 



(3) 



The self-energy £ has a normal (n) component at all 
T and and extra co component at T < T co . Its k- 
independence (locality) is the fundamental assumption 
of the dynamical mean field method. Off diagonal ~ r x 
components of E are not local in cordinate space, so do 
not appear in the DMFT. The two self energies are fixed 
by two coupled dynamical mean- field equations 14 for 
effective fields (a, b): 
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W 2 



(5) 



Here Gaa/bb are the two nonzero components of 
the local Greens function defined by G\ oc (i^>) — 
l/NJ2gG ( k, (recall J2k ek = 0- They are related 
to the effective fields {a n , b n } by Qaa (iw n ) — 9 1 * Za with 



Za= J drP({a n },r) and 



P(K},r)=, 



(6) 



(Ob and Zb analogous). P({a„},r) is the probability 
distribution of local distortions on the A-sublattice. 

For zero temperature we recover the previously dis- 
cussed commensurate charge order, because for Eq. [l| at 
T = Y^aa (k> ik>ri) = ~A and "Ebb (k^iu)^ = A. For 
T > we solve Eq. (Q) to @ using the procedure of Ref. 
|p"3f ; we retained 2 10 Matusbara frequencies and extended 
the frequency range analytically. We obtained conver- 
gence for T > 0.005 which is roughly T co for A = 0.64 
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Fig 1 Main panel: phase diagram of Hhol- Solid points: numer- 
ically calculated transition temperatures. Open points: 'modified 
BCS' approximation discussed below Dashed line: conven- 

tional BCS approximation. Shaded region: charge ordered; un- 
shaded: no order. Regime I: fermi liquid like normal state; regime 
III: 'polaron-like' normal state; regime II cross-over, as defined from 
P(r) and p(io) as discussed in text below. Left inset: T co vs A on 
logarithmic scale Right inset: gap to transition temperature ratio. 



Fig. 1 shows our calculated phase diagram. Where 
there is overlap (roughly A > 0.8) our results are in agree- 
ment with those of Ciuchi and de Pasquale H] if we scale 
T -> T/W and A -> \/{2W) with W = 2. From the 
left inset we see that the BCS approximation (dashed 
line) overestimates T c by a constant factor, but a differ- 
ent 'modified BCS' approximation (open circles) agrees 
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well with the small A numerical results. From the right 
inset of Fig. 1 we see that only for A — is the BCS 
value is obtained. 
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Fig. 2 Distribution of local distortions P (r) for three different 
coupling strengths A, for T = (vertical line), T m 0.93T CO (light 
lines) and T = 0.125 >> T co (heavy lines). For T < T co solid line: 
A sublattice; dashed line: B-sublattice. 

We next consider the phonon probability distribution 
P (r) and the local electron spectral function p (ui) — 
— — IuiGa/b shown in Figs 2 and 3 for several cou- 
plings and temperatures. 

A, = 1.44 A =2.40 A = 4.00 
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Fig. 3 Spectral function of electrons p (lj) for three different 
coupling strength A. for three different coupling strengths A, for 
T = (vertical line), T 0.93T CO (light lines) and T = 0.125 >> 
Tco (heavy lines). For T < T co solid line: A sublattice; dashed line: 
B-sublattice. 

For a weakly coupled fermi liquid (left panels, Figs 
2,3; region I of hgure I) P (r) is sharply peaked about 
its mean value r = and p (u>) takes approximately 
the non-interacting form (in our case, a semicircle). For 
a very strongly coupled system (right panels, Figs 2,3; 
region III, Fig. 1) one should think of the excitations 



as polarons, strongly bound electron-phonon complexes. 
In this situation P (r) displays two broad peaks (corre- 
sponding to distortions on occupied sites and antdistor- 
tions on unoccupied sites) and p (ui) has a two peak struc- 
ture corresponding to bound and empty sites. The cross- 
over between regions I and III occurs for couplings near 
A c = 37r/4 = 2.356 at which the weak coupling theory 
predicts a T = transition to a charge-disordered po- 
laron insulator [jl3| (this transition is preempted by the 
charge ordering transition). We have defined the left- 
hand boundary of the cross-over region as the locus of 
points at which P (r) loses its maximum at r = and 
the right-hand boundary as the locus of points at which 
p (u>) gains a minimum at uj = 0. The point A c sits ap- 
proximately in the center of region II. 

We now examine the weak coupling limit of T co ana- 
lytically. We first assume the equations have been solved 
for T > T co , yielding a mean field function c„. Just be- 
low the transition temperature the two effective fields a n 
and b n differ only slightly from each other or from the 
T > T co result, thus a n = c„ + e„ , b n = c„ - e n . 
Expanding the equations for the effective fields to linear 
order in e n we get the matrix equation for e„ 
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1 n €n 



Qn Gn 



Grri Gn 



(7) 

with Q n = j J dr P [j c 2g^ ■ T co is the temperature at 
which Eq. Mhas a solution. As written it is general; in the 

weak coupling limit G n = jr (l + A \ j^y + O (T 2 A 2 ) A 
withc° = (^) _1 = | (iu n + ^{iuj n ) 2 -W^j and 



A = A/ I H — — 



XW 
J 



E 



(c°; 



ee (A- 1 -A- 1 )" 1 (8) 



From the form of Gn we see that the expansion parameter 
is TA rather then A. The factor T is the classical ana- 
logue of the Migdal parameter ujd while the presence of 
A rather than A indicates that proximity to the polaronic 
instability strongly renormalizcs the coupling. Inserting 
G n into (0) and rearranging leads, after lengthy calcula- 
tion to the weak-coupling T co -equation 



1 = -A- 



1 



E 



G° 



(9) 



Eq. 



[3 ^ iuj n -2\T co Gl 
may be evaluated to logarithmic accuracy by 



replacing G by its uj = values ±2/W, yielding = 

7 7 W 

' , fr 



A (7 = 1.78107), showing explicitly the devia- 
tions from BCS. 

We now compare Eq |^ to the result obtained from the 
divergence of the charge order susceptibility defined by 
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E e_l 0>I Jdr((n t (r) n (0) ) - ( n (0) )f 



= ( 7? ) E Xnp Xnp ma Y be expressed in terms of an ir- 

reducible vertex function T which, within the dynamical 
mean- field approximation, is local fl5{|l6[| , thus Xnp sat- 
isfies Xnp = Xn P $np ~ Xn \ E 1 "nm Xmp with Xn = 



This equation 



may be expressed diagrammatically; in the weak coupling 
limit, the terms are shown in Fig. 4, which is exact up 
to terms of order (TA) 2 . 
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Fig. 4 Diagrammatic expansion of the one-particle self-energy 

(a) and of the irreducible vertex of the charge-order susceptibility 

(b) . Heavy lines: full local Green's function Q n . Dashed lines: free 
phonon propagator D^. All not shown but possible diagrams are 
smaller by additional factors of T. 



We find E r , 



AT and T I, 



-A + XSr, 



with 



A = A /(l + A 4 E Gn)- Note, A coincides with the pre- 
viously introduced one if we replace Q n — ► G n , which is 
correct at the order in T to which we work. Use of this 
f yields = Xo 1 ~ A with Xo = \ E l+XTxl '' ThiS 
diverges when Eq. || is fulfilled. 

We can now interpret Eq.(|j]) on more physical grounds. 
In the BCS limit, the irreducible vertex is given by A and 
X n by the convolution of two bare Green's function. The 
leading effect of increased coupling is a T-dependent self- 
energy, corresponding to inelastic scattering of electrons 
by phonons. Insertion of self-energies into x n increases 
the scale at which the logarithmic divergence is cut off, 
thereby reducing T co . This physics however does not op- 
erate at T = 0; thus it does not reduce the gap. In a more 
realistic model including quantum phonons some inelas- 
tic scattering would be present at u> > even at T = 0, 
renormalizing the gap slightly, but it would be weak in 
comparison to the thermally generated scattering so T co 
will be reduced more then A. Also in a more realistic 
quantum model the thermal scattering is negligible for 
temperatures lower than 0d/3 (Od is the Debye tem- 
perature); thus the effects we have described will be im- 
portant only if T co ^ 9d /3. In addition to the self-energy, 
there is a vertex correction, which acts just as the self- 
energy, to weaken the instability via inelastic scattering. 



The 'modified BCS' result ^| is shown as the dashed 
line in Fig. 1 and is seen to break down at the boundry 
of region I, when p (u>) and P (r) begin to acquire struc- 
ture indicating polaronic features. This shows that in the 
weak coupling regime, inelastic scattering is the key is- 
sue, and that 'pseudogap' formation not easily described 
in terms of traditional diagrammatics leads to changes in 
behavior. 

We finally discuss the relation of our results to the in- 
teresting proposal of McMillan jlo| concerning phonon 
entropy effects. As can be seen from Fig. 2, at interme- 
diate coupling thermal fluctuations of the phonon fields 
are crucial. If the fluctuations were negligible as in the 
BCS theory the widths of the peaks at T = 0.04 would 
be negligible in comparison to their separation. How- 
ever, as can also be seen the peaks either in the phonon 
distribution or the density of states do not shift appre- 
ciably until T becomes extremely close to T co . Thus the 
fluctuations, while strong, apparently do not feed back 
and reduce the mean value of the gap in the manner as- 
sumed by McMillan. However, the fluctuations do lead 
to a non-negligible probability of r « 0, i.e., of regions 
with a very small local gap. Electronic excitations in 
these regions give rise to the non-vanishing low energy 
density of states at the chemical potential, and as the 
regions become more probable, cause the destruction of 
the ordered phase. Comparison of the midlle and left 
panel of Figs. 2 and 3 shows that as the coupling is fur- 
ther increased, the physics crosses over smoothly to that 
expected in the very strong coupling limit. 
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